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ABSTRACT 

PDS 456 is a nearby (z=0.184), luminous (Z-boi ~ 10 47 ergs -1 ) type I quasar. 

A deep 190 ks Suzaku observation in February 2007 revealed the complex, broad 
band X-ray spectrum of PDS 456. The Suzaku spectrum exhibits highly statis- 
tically significant absorption features near 9 keV in the quasar rest-frame. We 
show that the most plausible origin of the absorption is from blue-shifted reso- 
nance (Is — 2 p) transitions of hydrogen-like iron (at 6.97 keV in the rest frame). 

This indicates that a highly ionized outflow may be present moving at near rela- 
tivistic velocities (0.26 — 0.31c). A possible hard X-ray excess is detected above 
15keV with HXD (at 99.8% confidence), which may arise from high column 
density gas (N\\ > 10 24 cm -2 ) partially covering the X-ray emission, or through 
strong Compton reflection. Here we propose that the iron K-shell absorption 
in PDS 456 is associated with a thick, possibly clumpy outflow, covering about 
20% of 4 7T steradian solid angle. The outflow is likely launched from the inner 
accretion disk, within 15-100 gravitational radii of the black hole. The kinetic 
power of the outflow may be similar to the bolometric luminosity of PDS 456. 

Such a powerful wind could have a significant effect on the co-evolution of the 
host galaxy and its supermassive black hole, through feedback. 
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1. Introduction 

Outflows are an important phenomenon in AGN and may play a key role in the co- 
evolution of the massive black hole and the host galaxy. Black holes grow by accretion 
and strong nuclear outflows can quench this process by effectively shutting off the supply 
of matter. At high redshifts, such quasar winds could have provided crucial feedback that 
controlled both the formation of stellar bulges and simultaneously self-regulated SMBH 
growth, leading ultimately to the observed M a relation for galaxies (Ferrarese k Merritt 
2000; Gebhardt et al. 2000). A wind, launched within a few gravitational radii of the black 
hole, can move into the host galaxy, pushing material outwards until, at a critical black 
hole mass, the ISM is effectively evacuated, quenching star formation and stopping further 
SMBH growth. Thus accretion disk outflows may help regulate the growth of the galactic 
bulge and the central black hole in galaxies (Silk k Rees 1998; Fabian 1999; King 2008; Di 
Matteo et al. 2005). 

Recently a number of high column density (N u ~ 10 2:5 - 10 24 cm~ 2 ), fast (> 0.1c) 
outflows have been claimed in several AGN (Chartas et al. 2002, 2003; Hasinger et al. 2002; 
Reeves et al. 2003; Pounds et al. 2003: Dadina et al. 2005; Gibson et al. 2005; Markowitz 
et al. 2006; Braito et al. 2007; Papadakis et al. 2007), through detections of blue-shifted 
K-shell absorption lines of iron, observed at rest-frame energies greater than 7 keV. These 
fast outflows are likely driven off the accretion disk by either radiation pressure (Proga et 
al. 2000; Proga k Kallman 2004; Sim et al. 2008) or magneto-rotational forces (Kato et al. 
2004), or both, a few gravitational radii from the black hole. In some cases the outflow rates 
derived can be huge, of several solar masses per year, equivalent to 1G 45 - 10 46 ergs _1 in 
kinetic power (Reeves et al. 2003; Pounds et al. 2003), matching the radiative luminosity of 
the AGN and can have a significant effect on the host galaxy evolution (King 2003). Such 
outflows are a likely consequence of near- Eddington accretion (King k Pounds 2003) and 
are key to understanding accretion at high Eddington rates. 

Here we present a deep 2007 Suzaku observation of the nearby (2 = 0.184, Torres et 
al. (1997)) and luminous (Ly ~ 10 47 ergs _1 , Simpson et al. (1999); Reeves et al. (2000)) 
type I quasar PDS 456. A previous 40 ks XMM-Newton observation revealed the presence of 
deep iron K-shell absorption above 7 keV, which may be attributed to a high column density 
outflow (Reeves et al. 2003), while PDS 456 also exhibits blue-shifted absorption in the UV 
(O'Brien et al. 2005). An analysis of the Suzaku XIS data reveals direct evidence for blue- 
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shifted absorption lines in the iron K-shell band (see Sections 3 and 4), which may arise from 
a near relativistic outflow, launched from the inner quasar accretion disk. We also report on 
a tentative detection of hard X-ray emission above 15keV in the Suzaku HXD observation 
of PDS 456. Values of H 0 ~ 70 kins 1 Mpc 1 , and Qa„ = 0-73 are assumed throughout and 
errors are quoted at 90% confidence (Ay 2 = 2.7), for 1 parameter of interest. 


2. Suzaku Observations of PDS 456 

PDS 456 was observed by Suzaku (Mitsuda et al. 2007) at the XIS nominal pointing 
position between 24 Feb to 1 Mar 2007 (sequence number 701056010), over a total duration 
of 370 ks. A summary of observations is shown in Table 1. Data were analysed from the 
X-ray Imaging Spectrometer XIS (Koyama et al. 2007) and the PIN diodes of the Hard 
X-ray Detector HXD/PIN (Takahashi et al. 2007) and processed using v2 of the Suzaku 
pipeline. Data were excluded within 436 seconds of passage through the South Atlantic 
Anomaly (SAA) and within Earth elevation angles or Bright Earth angles of < 5° and < 20° 
respectively. 


2.1. XIS analysis 

XIS data were selected in 3x3 and 5x5 editmodes using grades 0,2, 3, 4, 6, while 
hot and flickering pixels were removed using the SISCLEAN script. Spectra were extracted 
from within circular regions of 2.9' diameter, while background spectra were extracted from 
4 circles offset from the source and avoiding the chip corners containing the calibration 
sources. The response matrix (rmf) and ancillary response (arfs) files were created using 
the tasks xisrmfgen and xissimarfgen, respectively, the former accounting for the CCD 
charge injection and the latter for the hydrocarbon contamination on the optical blocking 
filter. Spectra from the two front illuminated XIS0 and XIS 3 chips were combined to 
create a single source spectrum (hereafter XIS FI), while data from the back illuminated 
XIS 1 chip were analysed separately. Data were included from 0.45 lOkeV for the XIS 
FI and 0.45 7 keV for the XISl chip, the latter being optimized for the soft X-ray band. 
The net background subtracted source count rates were 0.521 dt 0.002 ctss -1 for XIS FI 
and 0.413 ± 0.002 ctss" 1 for XISl, with a net exposure of 179.7ks. XIS background rates 
correspond to only 1.5% and 3.3% of the net source counts for XIS FI and XIS 1 respectively. 
The XIS spectra were subsequently binned to a minimum energy width corresponding to the 
HWHM XIS resolution (~ 65 eV at 6keV). Channels were additionally grouped to achieve 
a minimum S/N of 5 per bin over the background and y 2 minimization was used for all 


subsequent spectral fitting. 

2.2. HXD analysis 

The Suzaku HXD/PIN is a non- imaging instrument with a 34' square (FWHM) field 
of view. For weak hard X-ray sources (i.e. many AGN) where the source count rate is 
only a few percent of the background rate, it is crucial to correctly subtract the variable 
non X-ray background (i.e. the particle or detector background) to obtain an accurate 
estimate of the net source count rate. The HXD instrument team provide two Non X-ray 
Background (NXB) events files for this purpose; the background A model (also known as the 
“quick” background) and the background D model (or “tuned” background). The systematic 
uncertainty of the background D model is believed to be lower than for background- A, 
typically ±1.3% at the ler level for a net 20 ks exposure (Fukazawa et al. 2009). 

To determine which model provides the most reliable estimate of the background, the 
background A or D lightcurve was compared to the lightcurve obtained from the Earth 
occulted data (from Earth elevation angles ELV < —5), in the 15 50keV band. The Earth 
occulted data gives a representation of the actual NXB rate, as this neither includes a 
contribution from the source nor from the Cosmic X-ray background (CXB). The Earth 
occulted data were corrected for deadtime, using the pseudo events file. Figure 1 (top and 
middle panels) shows the comparison for each of the background models with the Earth rate. 
Specifically, we calculated the ratio (Earth - bgd [A or D]) / bgd [A or D] ; i.e. the difference 
between the Earth and background model rates normalized to the background model rate. 
The lightcurves are heavily binned (15 satellite orbits or 86.4 ks) in order to increase signal to 
noise. For the background A model, the model clearly underpredicts the Earth background 
rate, by +8.8 ± 1.0%. In contrast for background D, the model slightly overpredicts the 
Earth rate, by the ratio — 2.4 ± 1.0%. Thus subsequently we adopt the background D model 
as the likely more accurate background model, noting that for this observation it slightly 
overpredicts the Earth background rate, by 2.4%, which we correct for in the subsequent 
analysis. 

The background D NXB events file was then used in conjunction with the screened 
source events file (from rev2 of the pipeline processing) to create a common good time inter- 
val. The files were thus used to calculate lightcurves (binned into 5 orbital bins or 28.8 ks) for 
both the source (deadtime corrected) and for the detector (NXB) background; note a ler sys- 
tematic error of ±1.3% was included in the NXB rate. Figure 1 (lower panel) shows the PIN 
source lightcurve (after subtraction of the NXB) from 15 50 keV. The expected contribution 
from the diffuse (non- variable) Cosmic X-ray Background (CXB), using the spectral form 
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of Gruber et, al. (1999), was also calculated; in the 15 50 keV band this is expected to con- 
tribute 0.017 ctss -1 to the total count rate (or a 15-50 keV flux of 7.9 x 10~ 12 erg cm" 2 s -1 ). 
We include a ±10% uncertainty in the absolute normalization of the CXB. Note a more 
detailed discussion of the CXB and possible contamination in the HXD field of view is given 
in Appendix A. 

Thus the mean net source count rate, after subtraction of both the NXB and CXB 
components (i.e. total background = NXB ± CXB), is (1.5 ± 0.3) x 10~ 2 counts s _I , and 
is about 5% of the total background count rate (0.33 counts s -1 ). This corresponds to a 
source flux of (7.3 ± 2.4) x 10" 12 erg cm -2 s' 1 * * (including 1.3% systematic uncertainty on the 
background) from 15-50 keV, assuming a power-law photon index of T = 2. In comparison 
the source flux derived from using background A is higher, (1.2 ± 0.3) x 10~ n erg cm -2 s" 1 
from 15 — 50keV; hereafter we adopt the more conservative flux derived from background- 
D. This represents the first possible detection of PDS456 in the hard X-ray band above 
15keV, at the significance level of 7cr (statistical) or ~ 3cr (statistical±systematic) above the 
background and is below the expected sensitivity level of the Swift-BAT all sky survey, of 
about 10“ 11 erg cm -2 s” 1 (Tueller et al. 2008). 

A net (deadtime-corrected) source spectrum was subsequently extracted in the 15-50 
keV range (subtracting both the NXB and CXB) which yielded a net PIN exposure of 
164.8 ks and the same net count rate as above from 15-50 keV. The deadtime rate for the 
source events was 6.7%. The response file ae-hxd-pinxinome3-20080129.rsp provided by the 
instrument team was used in all the spectral fits below. The HXD/PIN spectrum was initially 
binned to a single bin, to represent a flux point between 15-50 keV, as described above. A 
finer binning of 2cr per bin above the background was subsequently used for the spectral fits 
in sections 4 and 5. Note that a cross-normalization factor of 1.16 is included between the 
HXD/PIN and XIS, as calculated for observations of the Crab. 1 


3. Detection of the Iron K-shell Absorber 

Initially we consider the X-ray spectrum of PDS 456 in the 0.45-10 keV band using the 

XIS alone. Galactic absorption of 2 x 10 21 cm -2 (Dickey k Lockman 1990; Kalberla et al. 

2005) was adopted, modeled with the “Tuebingen Boulder” absorption model (TBABS in 

XSPEC) using the cross-sections and abundances of Wilms et al. (2000). Figure 2 (upper 

panel) shows the residuals to the XIS spectra to a single power-law of T = 2.35 ± 0.01 with 

Galactic absorption. Evidently this is a poor fit (fit statistic, * 2 /dof = 1004.4/391); the 


0.45 lOkeV spectrum is concave, steepening towards lower energies, while there appears 
to be strong soft X-ray line emission observed near 0.75 0.80 keV and absorption at higher 
energies above 7.5 keV. 

We thus proceed to parameterize the XIS continuum with a broken power-law with 
r^ft = 2.44 ± 0.04, a break at E — 2.0 ± 0.2 keV and Fh^d = 2.24 ± 0.02 towards higher 
energies. Soft X-ray line emission is required, with one strong and broadened line at E = 
0.91 ± 0.01 keV (QSO rest frame) with an equivalent width of EW = 42 ± 10 eV, intrinsic 
width of a = 75 ± 15 eV (FWHM 57000 kins " 1 ) and a weak line at E = 1.16 ± 0.02 keV 
with EW = 6 ±3eV. Note a summary of the spectral fit parameters is shown in Table 2. 
The former line may be associated with either the NelX triplet (0.909 0.921 keV), L-shell 
(3d—2p) emission from Fexix (at 0.92 keV) or a blend of the two, the latter may be identified 
with FeXXiV L-shell (1.16 keV). 

The fit statistic is still poor (x 2 /dof = 529.3/382, null hypothesis probability of 8.3 x 
10" 7 * * ) and is significantly improved by the addition of two blue-shifted absorption lines ob- 
served at 7.67 ± 0.04 keV and 8.13 ± 0.07 keV (or at 9.09 ± 0.05 keV and 9.64 ± 0.08 keV in 
the quasar rest-frame), with A* 2 = 41.1 and 29.0 for each line. The line EWs are high at 
—137 ± 59 eV and — 135 ± 55 eV respectively. The intrinsic width of the absorption lines is 
loosely constrained to be a = 105^1 eV, thus the 90%> upper-limit on their FWHM velocity 
is < 12000 km s -1 (note the XIS resolution at 8keV corresponds to ~ 5000 km s ~ 1 FWHM 
at 8keV). Note the width of the lines are tied to have the same value as each other. 

The possibility that the absorption arises from a single, but broader line was also inves- 
tigated. We note that although this cannot be firmly ruled out, the resulting fit is slightly 
worse (Ax 2 = 12 for 1 parameter) than for the fit with 2 absorption lines, while the line 
width is substantially broadened (cr = 0.40 ± 0.12 keV or a FWHM of ~ 30000 kins' ‘). 

There is also evidence for broadened ionized iron Ko emission at 7.2 ± 0.2 keV in the 
quasar rest-frame (EW = 132 ± 52 eV, a = 0.35^;^ keV or ~ 30000 kms" 1 (FWHM), 
Ax 2 = 27.7). The final fit statistic is then x 2 /dof = 431.5/374. We will return to the origin 
of the iron K and L-shell emission in Section 5. 

We note that the absorption line detections are robust against the detector background 
of the XIS FI CCDs, which is < 2%, of the net source count rate, while the strongest 
emission line in the background spectrum (Ni I Ko at 7.47 keV) does not coincide with either 
of the absorption lines (note that unlike in XMM-Newton EPIC-pn spectra, no Cu Ko line 
is present in Suzaku XIS background spectra at 8keV). The comparison between the net 
source spectrum and background spectrum is plotted in Figure 3. 


See ftp: / /legacy.gsfc.nasa.gov/siuaku/doc/xrt /suzakumemo-2008-06.pdf 
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3.1. Statistical Significance of the Absorption Lines 

The significance of the iron K absorption lines was estimated using multi-trial Monte 
Carlo simulations. The method used was identical to that described in Porquet et al. (2004) 
and Markowitz et ah (2006). Specifically we generated 1000 random XIS-FI spectra with an 
initial input continuum, exposure and background as per the actual PDS 456 data, but under 
the null hypothesis of no absorption lines. The spectra were stepped in energy increments 
corresponding to the HWHM XIS energy resolution over the 3 lOkcV energy band, i.e. where 
any red or blue- shifted Fe absorption lines may be reasonably expected to be found, in order 
to map the distribution of A\ 2 for the trials. In the actual data, two absorption lines were 
found, yielding improvements to the fit statistic of Ay 2 = 41.1 and 29.0 respectively. In the 
simulated spectra, these values correspond to 3 and 7 out of the 1000 spectra respectively. 
Thus the null hypothesis probability that each of the individual lines is due to random 
statistical fluctuation is 3 x 10~* and 7 x 10“ 1 respectively. Furthermore the null probability 
that both lines are random is the product of these probabilities, i.e. 2 x 10 -5 . 

We further note that the 9 keV rest-frame absorption is also independently detected in 
both XIS FI CCDs, XIS 0 and XIS 3, but not the back-illuminated XIS1 which has little 
sensitivity above 7 keV. The likelihood that the absorption is due to random chance is further 
diminished when the 2001 XMM-Newton data are considered. Figure 2 (lower panel) also 
plots the 2001 EPIC-pn data overlaid upon the Suzaku XIS spectrum. Absorption is clearly 
present in the XMM-Newton data (with a trough apparent between 9-10 keV, quasar rest 
frame), although the exact absorption profile may have varied. Indeed the possibility of 
deep iron K absorption in PDS 456 was first noted in Reeves et al. (2000), on the basis of 
low resolution RXTE PCA data. Thus the coincidence of absorption in two independent 
observations (and telescopes) at very similar energies would appear to confirm the detection 
of the blue-shifted absorption in PDS 456. 


4. The Origin of the Iron band Absorption 

The possible origin of the iron K band absorption is now discussed. 


4.1. Identification of the Absorption Lines 

As shown above, the absorption lines detected in the XIS appear to be highly blue-shifted 
with respect to the iron K-shell band, with QSO rest frame energies of 9.09 ± 0.05 keV and 
9.64 ± 0.08 keV. Figure 2 (lower panel) shows the XIS FI residuals in the iron K band to the 
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above continuum model in the quasar frame. Thus if either line is associated with the strong 
Fexxvi Lyman- o transition at 6.97 keV, then this implies relativistic velocity shifts of v on% \ = 
—0.26 ± 0.01c (—78000 ± 2000 kms -1 ) and u out2 = —0.31 ± 0.01c (—93000 ± 3000 kms -1 ) 
respectively for each line, while if the absorption is associated with lower ionization Is — 2p 
transitions, e.g. FeXXV at 6.70 keV, the blue-shift required will be correspondingly higher. 
Moreover an identification without a large blue-shift is extremely unlikely; there are no 
strong line transitions in the 9.T 9.6 keV range from abundant elements (e.g. ZnXXX ls-2p 
at 9.3 keV or Nixxvm Is - 3 p at 9.58 keV would be undetectable barring extraordinarily 
strange abundances). Bound free absorption by Fexxvi at 9.27 keV may appear to be 
plausible, however in a self-consistent photoionization model (see below), a strong (> 100 eV 
EW) Fexxvi Ka absorption line (at a rest-energy of 6.97 keV, 5,S9keV observed) would 
accompany a detectable (i.e. r > 0.1) edge-like structure at 9.3 keV; such a line is not 
present in the spectrum to a limit of < 20 eV in EW. We discuss this possibility further 
below. 

Neither is it likely that the lines arise from local ( z = 0) ionized matter, for example 
associated with our galaxy or the local IGM; e.g. McKernan et al. (2004, 2005), but also see 
Reeves et al. (2008). In this scenario the line at 7.67 keV (observed frame at z = 0) could 
only arise from K/3 transitions of Fexxm or Fexxiv from 7.6-7. 8 keV, while the 8.13keV 
(observed) line could be associated with either the Kfi line from Fexxvi (8.25 keV) or the 
Ly-a line of Nixxvm (8.10keV). However the K/3 lines of ionized iron (i.e. from Fexvm 
and above) can be excluded, as the correspondingly stronger Ko transitions between 6.5 
7.0 keV are not observed, while any Ni K-shell absorption will be extremely weak due to its 
low abundance. 

Furthermore any absorption from the local IGM would be of rather low turbulence 
velocity. For instance if a component of the local IGM or local hot bubble gas has a temper- 
ature of T — 10 7 K, then the thermal width of the lines would only be a = 40 kms -1 for iron 
atoms. Indeed any absorption from this gas with such a low thermal width would produce 
very weak absorption lines, with equivalent widths of the order a few eV, undetectable with 
the current non-calorimeter detectors. Thus local ( z = 0) matter can be ruled out as a cause 
of the 9keV absorption. Instead the absorption must be associated with PDS 456, possibly 
from highly blue-shifted absorption from an outflowing wind, as is discussed next. 

4.2. Photoionization Modeling of the Iron K Absorption 

We further model the Fe K absorption using a grid of XSTAR (v2.11n9) photoionization 
models. This uses the latest Fe K-shell treatment of Kallman et al. (2004) and solar abun- 
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dances (Grevesse &: Sauval 1998), while an illuminating continuum from 1 1000 Rydbergs 
of T = 2.2 is assumed. A large turbulence of > 1000 kins -1 is required in order for the 
Fe K absorption lines to be un-saturated and to be of high (~ 100 eV) equivalent width; 
such a velocity may be associated with an inner accretion disk wind. Here we adopt a high 
turbulence velocity of er tU rb = 10000 km s -1 , consistent with the limits on the absorption line 
width, noting that we also tested the models (A C) below using a lower turbulence grid of 
3000 km s -1 , producing almost identical results. However smaller turbulences < 1000 km s -1 
fail to reproduce the observed absorption line equivalent widths due to saturation. The XIS 
FI only is used in the spectral fits, as this is most sensitive towards the absorption above 
7keV. 

The continuum is parameterized by a broken powerlaw absorbed by a Galactic line of 
sight column of gas over the 0.5-10 keV band, as discussed previously. The two soft X-ray 
Gaussian lines that are apparent near 1 keV, are also included; hereafter we use this as the 
baseline continuum model for comparison. Although this provides a good description of the 
broad-band X-ray continuum in the XIS FI, the data show clear residuals from 6 -9keV in 
the iron K band (e.g. Figure 2). The fit statistic of * 2 /dof = 298.8/210 is formally rejected 
with a null hypothesis probability of 5.4 x 10~ 5 . 

The required iron K absorption is first modeled with a single layer of photoionized gas 
fully covering the X-ray source, fitted via a grid of absorption models generated by Xstar 
and is hereafter referred to as Model A. We can express the model as:- 

F(E) = (BPL + GA em ) X XSTAR(£, N H , tw) x (1) 

where BPL is the broken powerlaw continuum, GAe m represents the Gaussian emission lines, 
XSTAR is the single layer of absorbing gas modeled via Xstar (with JV H , £ and v oat as variable 
parameters), e ~ a(E)NHG * 1 is the neutral Galactic line of sight photoelectric absorption and 
F(E) is the emergent energy spectrum. This produces a good fit to the data, (y 2 /dof = 
224.0/204; null probability 0.151), a blend of the FeXXV and Fexxvi ls-2p resonance lines 
models the absorption line profile; the Xstar absorber parameters are summarized in Table 3. 
A large outflow velocity of v out = —0.30 ± 0.02c is required in order to reproduce the energy 
of the absorption lines (near 9keV in the rest frame). Furthermore a high column density 
(Ah = 2.2^} o x lO^cm -2 ) and ionization parameter (log£ = 4.1 ± 0.2 2 ) are also required 
to produce strong highly ionized K-shell absorption lines from iron. The iron K-emission is 
included as a simple broad Gaussian, noting that we will consider more physical models for 
the iron K/L shell emission in Section 5 and 6. 


2 The units of £ are erg cm s K 
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However one Xstar zone with a single outflow velocity does not quite accurately model 
the absorption profile shown in Figure 4a, as the spacing between FeXXV and Fexxvi 
Is — 2p transitions (6.70 keV and 6.97 keV) is smaller than the observed separation of the 
absorption lines at 9.08 keV and 9.62 keV (rest-frame) respectively. Thus a better fit is 
subsequently obtained using two zones of absorbing matter, with different outflow velocities 
but assuming identical columns and ionizations. This model is referred to as model B, and 
is shown in Figure 4b. A column density of N H = (1.3 ± 0.5) x 10 24 cm 2 is required, while 
a high ionization parameter of log£ = 4.91J;? implies that iron is predominantly in a H-like 
state. The outflow velocities obtained are v oiitl = -0.26 ± 0.02c and v out2 = - 0.31 ± 0.02c, 
while the fit statistic is good (* 2 /dof = 207.3/203, null probability 0.4). Note that the 
presence of substantial columns of lower ionization gas (with log £ < 3) fully covering the 
AGN are excluded, as this would produce considerable absorption from L-shell iron, as well 
as K-shell lines from Mg, Si and S below 3 keV, which are not observed in the Suzaku 
XIS spectrum. Low ionization gas, partially covering the emission below 10 keV during the 
Suzaku observation, can be excluded due to the lack of convex curvature in the XIS spectrum; 
however this is not the case for the 2001 XMM-Newton spectrum of PDS456 discussed in 
Section 5.3. 

We also allow both the iron emission and absorption to have the same self consistent 
velocity broadening, via convolution with a simple Gaussian broadening function, hereafter 
referred to as model C. This might be the case kinematically if the line emission is integrated 
over a smooth wide angle outflow or if the absorption is viewed over multiple lines of sight, e.g. 
if the absorber is extended w.r.t the X-ray source, or if photons are scattered throughout the 
flow, as might be expected for the high columns observed here; This model can be expressed 
as:- 

F(E) = [(BPL + GAen.) x XSTAR (£, A„, t; out )j x GA conv x r r^ E)N ^ (2) 

where GAconv is the Gaussian broadening convolution. For simplicity a single bulk outflow 
velocity is assumed for the absorption, while the iron emission / absorption is assumed to be 
intrinsically narrow (prior to broadening), except for the turbulence velocity noted above. 
The resulting model is shown in Figure 4c, which has a velocity broadening of 26000 ± 
7000 km s -1 (FWHM) and a net outflow velocity of —0.29 ± 0.02c. While this represents a 
statistically adequate fit (* 2 /dof = 227.9/205, null probability 0.131), the model does not fit 
the profile of the absorption feature at 8 keV as well (c.f. model B), which appears narrower 
than the velocity broadening derived here. This might suggest that a homogeneous outflow 
is too simplistic and that significant clumping in the absorber inav occur; this is discussed 
further in Section 6. 

Alternatively we attempt to model the iron K absorption with no net outflow velocity, 
in the quasar rest, frame at z = 0.184. This is only plausible for very low values of the 
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turbulence velocity (e g. o tm u = 100 km s -1 or less), so that the 9 keV rest-frame absorption 
is modeled by deep bound-free edges of FeXXV and FcXXVl, yet the equivalent widths of 
the corresponding I s — 2 p absorption lines are small, as no absorption is observed in the 
spectrum at 6.70 or 6.97 keV (5.66, 5.89 keV observed frame). Thus the high turbulence, 
high outflow velocity absorber from models A or B above is replaced by a low turbulence 
(cTtnrb — 100 km s ' 1 ) absorber, with no net outflow velocity. The resulting model is listed 
as model D in Table 3. A very high column density of N H = 2.3 + q 2 x 10 24 cm -2 and high 
ionization of log£ = 4.3l[| 2 is required in order to produce a deep absorption edge near 9 
keV. Nonetheless the fit obtained with the low turbulence absorber is significantly worse 
than for the high velocity absorber (\' 2 /dof = 254.0/208, null probability 1.6 x 10~ 2 ) and the 
model only crudely fits the profile of the absorption, as shown in Figure 5. In addition even 
for the low turbulence of 100 kins” 1 , the model still predicts significant resonance Is — 2 p 
absorption between 5.6-6.0keV in the observed frame, which is not present in the data. 

Thus the most plausible origin of the absorption appears to be from highly blue-shifted 
iron K -shell lines of FcXXVl or FeXXV associated with a mildly relativistic (0.25-0.3c) out- 
flow. Hereafter we adopt model B is the best parameterization of the iron K band absorption. 

5. Modeling the Broad-Band Spectrum of PDS 456 

PDS 456 appears to be weakly detected in the hard X-ray band, as described in Section 2, 
with an integrated 15 5QkeV band flux of 7.3±2.4xl0 -12 ergcm~ 2 s~ 1 . This flux is somewhat 
higher than the flux measured by XIS (3.8 ±0.1 x 10~ l2 ergcm" 2 s _1 ) in the 2-10 keV band. 
Indeed if we extrapolate the best fit XIS broken powerlaw continuum model in the previous 
section to higher energies, the predicted 15-50 keV flux is only 2.1 ±0.2 x 10~ 12 erg cm -2 s -1 . 
This possible factor of ~ 3 increase in flux above 10 keV is illustrated in Figure 6, which shows 
the extrapolation of the continuum fitted to XIS out to higher energies. The absolute level 
of the HXD flux is dependent on the background systematics of the detector, as discussed in 
Section 2, which we account for as an additional uncertainty on the HXD/PIN background 
level. After inclusion of this systematic error, the excess in HXD flux in Figure 4 appears 
significant at a moderate level of 99.8% confidence (corresponding to Ay 2 = 13), or about 
3cr . Further Suzaku HXD observations will hopefully confirm the validity of the hard X-ray 
detection of PDS 456. 


5.1. Reflection Dominated Models 

For now, we proceed by taking the HXD/PIN flux measurement at face value, noting 
that the fits below include the systematic error on the HXD/PIN background. Compton 
down-scattering or “reflection” off optically-thick matter (Light man k White 1988; George 
k Fabian 1991), e.g. the disk or torus, could in principle account for the hard X-ray excess. 
We model this with Compton reflection off partially ionized material, using the REFLION 
model (Ross et al. 1999; Ross k Fabian 2005). 

The emergent reflected spectrum is calculated for an optically-thick atmosphere (such as 
the surface of an accretion disk) of constant density illuminated by radiation w'ith a power- 
law spectrum, with a high-energy exponential cutoff with an e-folding energy fixed at 300 
keV, i.e. varying as e - £ / 300keV . The reflected spectrum is calculated over the range leV 
to 1 MeV. Non-LTE calculations provide temperature and ionization structures for the gas 
that are consistent with the local radiation fields. In addition to fully- ionized species, the 
following ions are included in the calculations: Cm- VI, Nm-vil, O III- VIII, Nelll-X, Mg III- 
XII, Siiv-Xiv, Siv-XVl, and FeVl-XXVl. Thus the model self consistently computes the iron 
line (K or L-shell) emission, as well as from other elements. 

The REFLION model is also convolved with relativistic blurring from a disk around 
around a black hole (Fabian et al. 1989; Laor 1991), to describe the reflected emission 
expected from the innermost accretion disk. A disk emissivity falling as R~ 3 is assumed, 
along with an outer disk radius of 400i?g, where R g is the gravitational radius. The model 
can be expressed as:- 

F(E) = [CPL ± (REF x KDBLUR)] x XSTAR(C N H , v out ) x c“ <7(E)Nf, G ^ (3) 

where CPL is the illuminating cut-off powerlaw continuum, REF is the Compton reflected 
emission convolved with relativistic blurring (KDBLUR). Absorption from the putative high 
velocity outflow (denoted as XSTAR), is also included (modeled by an Xstar grid as per 
model B as described in Section 4.2) in order to account for the iron K-shell absorption at 
9keV. 

An adequate fit (y 2 /dof = 248.5/212, null probability 4.3 x 10~ 2 ) is obtained, with 
an inclination angle of 56 ± 8°, an inner radius of < 20/? g , an ionization parameter of 
£ = 190 ± 95 and a power-law photon index of T = 2.32 ± 0.03. The reflection fraction is 
then R = Q/2n = 1.3, where R = 1 corresponds to an infinite slab effectively subtending 
27 r steradian to the X-ray source. 3 A high Fe abundance of x4 Solar is found, driven by 


3 Note that we do not calculate the error on R, as in the REFLION model the reflected ratio between 
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the strong iron L-shell emission near 1 keV, although this value is not well constrained. 
Overall the fit is good below 10 keV, however the model fails to reproduce all of the excess 
in the HXD above 10 keV, see Figure 7, although the precise HXD/PIN flux is subject to 
systematic uncertainty as discussed above. Thus a contribution from a reflection component 
towards the X-ray emission from PDS456 may be present, but does not appear to account 
for all of the hard X-ray emission observed in the HXD. 


5.2. Absorption Models 

An alternative scenario is that the majority of the intrinsic continuum flux from PDS 
456 is in fact absorbed below 10 keV. therefore the heavily absorbed continuum only emerges 
above lOkeV and accounts for the observed hard X-ray excess. This is similar to the scenario 
in Seyfert 2s, where a heavily absorbed high energy continuum is often seen (Risaliti 2002), 
except in PDS 456 the absorber must only partially cover the continuum X-ray emission in 
order for some direct flux below lOkeV to leak through, producing the observed continuum 
flux in the XIS. A similar, Compton-thick, but partially covering absorber has also recently 
been detected in the Suzaku HXD observation of the luminous AGN (at z = 0.104), 1H 0419- 
577 (Turner et al. 2009). Any partial covering absorption must be located close to the central 
AGN, i.e. much closer than any pc scale absorber, in order to partially cover the compact 
X-ray emission. 

In order to model this absorption we include a second layer of photoionized gas via the 
same Xstar grid as described in Section 4.2 (with solar abundances and cr tU rb = 10000 kms" 1 ), 
in addition to the highly ionized absorption that is responsible for the iron K-shell lines at 
9 keV. Thus this second layer of photoionized gas absorbs only a fraction / of the intrinsic 
cut-off power-law continuum, while a fraction 1 — / is not absorbed. Hence the model can 
be expressed as:- 

F(E) = [(CPL x / x PC) 4- (CPL x (1 -/JJ+GAo,,,] xXSTAR(f, iV H , iw) x e _ ' ,lE|A,HG * 1 (4) 

where CPL is the cut-off powerlaw continuum with an e-folding energy of 300 keV, PC is 
the layer of absorbing gas that partially covers the X-ray source with a covering fraction /, 
while XSTAR represents the highly ionized outflowing gas that is responsible for the 9 keV 
iron K-shell absorption, which we assume to fully cover all of the X-ray emission. In this 
model the iron K and L-shell emission is included as Gaussian line components represented 
by GAem 


the reflected emission and the power-law continuum is calculated over the range 1 eV to 1 MeV, with large 
uncertainties in the extrapolation. 


This model provides a statistically better fit (y 2 /dof 217.9/209, null probability 
0.32), with a photon index of F = 2.27 ± 0.03. A large column density is required to model 
the heavily absorbed component above lOkeV, with N n - (3.5 ± 1.5) x 10 24 cm“ 2 and an 
ionization parameter of log£ = 2.5 ± 0.4 erg cm s" 1 ; this lower ionization parameter value 
(log£ < 3.0) is required for the resulting continuum to be largely opaque below lOkeV, with 
most of the opacity arising from partially ionized iron (i.e. FeXXIv and below). Likewise this 
sets the lower limit on the column density (jV H > 2 x 10 24 cm' 2 ). An outflow velocity for the 
partial covering absorption zone is not formally required, however it is not well constrained 
with a lower limit of < 0.2c, less than for the highly ionized zone. The (line-of-sight) covering 
fraction of the absorber is f mv = 0.724:0.12, i.e. ~ 30% is unabsorbed and this is the emission 
viewed below lOkeV. 

Note that one drawback of the partially covering model discussed above, compared to 
the reflection model, is that the iron K and L-shell emission is not computed, instead we have 
simply parameterized the emission by Gaussi&ns as discussed in Sections 3 and 4. However 
in the case of a Compton thick wind (i.e. with N\\ > 10 24 cm 2 ), significant emission via 
Compton reflection off the wind would be expected and could account for observed emission 
lines as well as the iron K band absorption. We discuss such a model further in Section 6.3. 


5.3. Comparison with Previous Observations 

Variations in the absorption towards PDS 456 may also explain the drastic long term 
spectral variability of the AGN. Figure 8 shows the uF v spectral comparison between the Feb 
2007 Suzaku observation with those by RXTE (1998/2001), ASCA (1998), XMM-Newton 
(2001) and Chandra (2003); see Table 1 for a list of observations. 4 Below 10 keV, PDS 456 
shows remarkable spectral variability. RXTE 1998 and XMM 2001 caught the AGN at a 
relatively high flux, as previously discussed. Reeves et al. (2003) derived a covering fraction 
of 60% for the absorber in the 2001 XMM-Newton observation, with a column density of 
5.71^5 x 10 23 cm" 2 , an ionization parameter of log£ = 2.5 jfc 0.3 and an outflow velocity 
of ~ 50000 km s" 1 , although the latter value was not well constrained in the shorter 40ks 
observation. The low resolution RXTE/PCA spectra, from both 1998 and 2001, also show 
a deep absorption trough at 8keV in the observed frame, with r ~ 1 implying a column 
density for the absorber of ~ 10 24 cm -2 (Reeves et al. 2000). 

We re-analysed the XMM-Newton EPIC-pn spectrum from 2001, using the same Xstar 


4 For clarity we show only the 1998 RXTE observation in Figure 8, the 2001 observation caught the source 
at a similar flux level. 
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grids as used in this current paper (with 10000 km s“ 1 turbulence). The current grids have the 
advantage in that they include up to date opacities from partially ionized iron (Kallman et 
al. 2004); e.g. Is 2 p transitions from Foxvni-xxiv and high order transitions from Fexvii 
and lower. A model of the same form as the partial coverer in Section 5.2 was adopted. 
As per Suzaku XIS. a high ionization outflowing zone can fit the observed dip at 8keV; 
with log £ — 4.5 ± 0.5, N = 5 ± 2 x 10 23 cm~ 2 and c out ~ —0.30c. Interestingly partially 
covering absorption is also required in the XMM-Newton data to model the substantial 
continuum curvature present in the data from 1-10 keV (see Figure 8, green points), however 
the column density is lower ( A r n ~ 1 x 10 23 cm -2 ) compared to the partial coverer in the 
Suzaku observation (An > 10 24 cm~ 2 ). The partially covering gas is low ionization (logf < 
2.1) and does not produce strong Is — 2 p lines from Fexvm-xxiv, as the L-shell is filled, 
while its outflow velocity is not constrained by the pn data. The covering fraction of this 
absorber is measured to be 80%, very similar to the Suzaku value. 

In contrast to this, the Chandra/HETG observation of PDS 456 in 2003 (and also 
ASCA in 1998) showed a very hard (T = 1.3) low flux spectrum (i.e. XMM-2001, F0.5-10 = 
1.1 x 10" 11 erg cm" 2 s _1 ; Chandra 2003, Fq.s-io = 4.0 x 10~ 12 erg cm -2 s -1 ). Changes in the 
covering and/or column of the (lower ionization) absorbing matter could possibly reproduce 
this spectral variability. For instance if the column density of the partial covering absorber 
was lower during the 2001 XMM-Newton observation, compared to the Suzaku observation 
(e.g. N\{ ~ 1G 23 cm 2 vs > 10 24 cm ' 2 ), this would allow more of the direct continuum to be 
observed below 10 keV, resulting in a higher flux during the XMM-Newton observation and 
accounting for the convex spectral shape below 10 keV. The lowest flux Chandra observation 
could also contain a greater contribution from the reflected emission, however any such 
component would have to be broadened as no narrow iron K-emission is observed. The 
long-term spectral variability of PDS 456, e.g. in the context of a variable partial covering 
absorber, will be investigated in more detail in a subsequent paper (Behar et al. 2009, in 
preparation) . 


6. Discussion 

6.1. The Intrinsic Luminosity of PDS 456 

One possible consequence of the high column density absorption towards PDS 456 is 
that the intrinsic X-ray luminosity may be higher than inferred below lOkeV; the observed 
luminosity from 2 lOkoV is only 4 x 10 44 ergs -1 . In contrast the bolometric luminosity 
is substantially higher. From the optical spectrum of PDS 456 (Simpson et al. 1999) and 
correcting for the cosmology used in this paper, i^L^( 5100A) = 2 x 10 46 ergs -1 ; thereby 


assuming a typical relation between the 5100.4 luminosity and the bolometric luminosity 
(L bol ) of L bo i ~ 9 x i/L„( 5100A) (Kaspi et al. 2000), an estimate of L bo i = 1.8 x 10 47 ergs _1 
for PDS 456 is derived. As a consistency check, we calculated the integrated infra-red to UV 
luminosity (L IR _ UV ) based on our previous spectra in the IR (1 — lO^rn), optical (3000 — 
SOOOA) and UV (1200 - 3000A) bands (Simpson et al. 1999; O’Brien et al. 2005). The de- 
reddened Lir_uv is 1.0 x 10 47 ergs _1 , which defines a lower-limit on L bo i, as the EUV band 
emission is not included. Thus hereafter we adopt Lt K ,i = 2 x 10 47 ergs 1 as the bolometric 
luminosity of PDS 456. 

Hence the observed 2-10 keV luminosity measured by the XIS is only 0.2% of bolometric, 
meaning PDS 456 is X-ray faint compared to most quasars, where 3 5% of bolometric may 
be expected for a typical AGN SED (Elvis et al. 1994). This then suggests that part of 
the intrinsic X-ray emission may be hidden or absorbed, as discussed in the partial covering 
model, where 70% of the emission may be absorbed below 10 keV by a near Compton- thick 
layer of gas. Correcting for this putative absorption below lOkeV, the 2-10 keV luminosity is 
then 1.2 ±0.3 x 10 45 ergs -1 . Furthermore at the high column densities observed in PDS 456, 
the observed continuum X-ray flux may also be suppressed via by a factor e~ T , where r = 
N H cr T and a T — 6.65 x 10~ 25 cm 2 is the Thomson cross-section. Thus for a column of 
jV H = 2 x 10 24 cm~ 2 , derived from the Xstar fit to the absorber, then r = 1.3 and thus 
the X-ray luminosity corrected for scattering and absorption is ±2-10 = 4 x 10 45 ergs -1 , 
about 2% of the bolometric luminosity of PDS 456. This value is consistent with the typical 
~ 3% ratio between L 2 -io and L bo i measured for typical quasars (Elvis et al. 1994). This 
makes PDS 456 one of the most luminous known nearby quasars in the X-ray band similar to 
3C 273, which is also consistent with its large bolometric luminosity, a factor of xl.7 higher 
than in 3C273 in the optical band (Simpson et al. 1999). 


6.2. An Outflow Model for PDS 456 

Given the high ionization state (and possible partial covering) of the absorber in PDS 456, 
it is more likely to reside close to the X-ray emission region and not at parsec scales com- 
mensurate with the molecular torus, as predicted by AGN Unified schemes (Antonucci 1993) 
and as suggested for some Seyfert outflows (Behar et al. 2003; Blustin et al. 2005). 

Thus it is not likely to be the same matter responsible for the X-ray absorption towards 
Seyfert 2 galaxies, as PDS 456 is a classic luminous type I quasar with an unobscured view 
of the optical/UV BLR (Simpson et al. 1999; O’Brien et al. 2005). The X-ray absorber in 
PDS 456 is probably too highly ionized to cause significant absorption of the UV emission, 
although an outflowing absorption component to the Ly a line profile appears in the HST- 
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STIS spectrum ( v = 14000 — 24000 km s -1 ), as well as highly blue-shifted ( v — 5000 km s" 1 ) 
Civ emission (O’Brien et al. 2005). Thus one possibility is that the X-ray absorber is 
both part of a massive outflow, with density variations responsible for different gas layers. 
Dense (71 ~ ID 10 cm -3 ) clumps within the outflow can account for the possible Compton-thick 
absorption above 15 keV, which partially covers the X-ray emission, while the less dense (and 
highly ionized) fully covering gas is responsible for the strongly blue-shifted Fe K absorption 
lines. Indeed the high ionization gas may well shield the lower ionization matter. 

We first consider the case of a homogeneous radial outflow from PDS 456, which could 
exist as a spherical flow, but more likely as a bi-conical wind (Sulentic et al. 1998; Elvis 
2000). For simplicity, we assume that the outflow velocity is approximately constant on the 
compact scales observed here, although initial acceleration must occur close to the launch 
radius (i.e. a radiatively driven wind), with deceleration occurring at larger radii. From 
simple conservation of mass, the outflow rate is:- 

Mout = 4TrbnR 2 m p v 0{lt (5) 

(where nR 2 = L\ on/0- Here 6 is a geometrical factor (i.e. b = 1 for a full covering, 
homogeneous spherical outflow). 

In PDS 456 the ionizing luminosity, defined by Xstar from 1 — 1000 Rydberg with an 
input continuum of T = 2.2 is L io ti = 3 x 10 4o ergs~ 1 , correcting for Galactic absorption. 
This luminosity may indeed by higher by a factor of 3, e.g. if 70% of the emission below 
lOkeV is absorbed, as in the partial covering scenario. However in the calculations below 
we adopt L mn = 3 x 10 45 ergs _1 as a robust lower limit to this luminosity. 

As discussed in Section 4.2, the high ionization matter appears to be outflowing with 
Com = —0.26c, with an ionization parameter of log£ = 4.9; e.g. see model C, Table 3. Thus 
the outflow rate in PDS 456 is of the order M out = 6 x 10 27 6gs _1 or 1006 M G yr _1 . Even for 
a rather conservative value of 6 ~ 0.1, then M uut ~ 10 M©yr _1 . The outflow kinetic power 
is then simply E out ~ A/ out e 2 ut /2, which is then 2 x 10 47 6ergs _1 . Thus the kinetic output 
of the wind may be an appreciable fraction of the bolometric luminosity for PDS 456 for a 
reasonable value of b. Note that the kinetic power is very sensitive to n out , effectively varying 
as ?> 3 ut , but note that we adopt the lowest value of the velocity range v out = 0.26 - 0.31c 
derived in Table 3. 

There is no direct (e.g. reverberation) mass estimate for the black hole in PDS 456, 
however we can estimate its likely value from known scaling relations, derived from rever- 
beration methods, between the AGN BLR virial radius and black hole mass (Kaspi et al. 
2000; McLure k Jarvis 2002). If we adopt the relation derived in McLure k Jarvis (2002):- 

M bh = 4.74(AL 5i ^/10 44 ergs- 1 ) O(51±o lo [FWHM( J H/?)] 2 M© (6) 


where AZ^^ = 2 x 10 46 ergs 1 and FWHM (Hfi) = 3974 ± 764 kins 1 (Torres et al. 1997; 
Simpson et al. 1999), then the black mass is logA/ BH = 9.3 ± 0.4 M 0 for PDS 456. Alterna- 
tively, using the equivalent relation in Kaspi et al. (2000) gives logM B n = 9.5 ± 0.3 M . 
The Eddington-limited luminosity for an accreting black hole of mass A/ H h is L md = 
4:7rGMm v c/(j T , while the Eddington accretion rate is simply Mem = Law /(r;c 2 ), where 
77 is the efficiency of converting rest-mass into energy. For PDS 456 if Af BH = 2 x 10 9 M©, 
then L md = 2 x 10 47 ergs -1 and Mem — 50 M 0 yr -1 for 7/ = 0.06 (tl re maximum efficiency 
for accretion onto a Schwarzschild black hole). Thus the mass outflow rate is likely to be a 
substantial fraction of the total accretion rate. 


The escape radius of the outflow is simply /?,*,. ~ c 2 /v 2 R tJ > 15 R g (where Z? g = 
GMbe/c 2 is the gravitational radius), equivalent to R^ c = 5 x 10 15 cm for A/ HH = 2 x 10 9 M©. 
For the case of a homogeneous radial outflow, the density varies as n x R~ 2 and thus the 
ionization parameter £ = L/nR 2 is largely independent of the outflow radius. Hence the col- 
umn density viewing through the outflow line of sight down to a characteristic wind radius 
Z^wtiid is ; - 


Nh= [ n(R)dR = f (L^/SR^dR = L, 


( 7 ) 


Therefore for PDS 456, as jV h = 1.5 x 10 24 cm" 2 and log£ ■= 4.9, then i? wind = 3 x 10 l6 cm ~ 
100/? g . Note R wlud is within the expected UV/BLR emission radius for PDS 456; e.g. for 
FWHM(Ly Q ) = 12000 km s' 1 (O’Brien et al. 2005), then /?uv ~ lQ0Q /? g . Note that as the 
highly ionized matter is located close to the black hole in PDS 456, this likely excludes the 
possibility of low turbulence (and low outflow velocity) gas (model D, Section 4.2) contribut- 
ing towards the iron K band absorption. 


It is also possible that the wind is not homogeneous but is instead clumpy. Indeed X-ray 
variability is observed in PDS 456 on rapid (20 — 30 ks) timescales, e.g. as seen previously 
in the XMM-Newton, Beppo-SAX or RXTE observations (Reeves et al 2000, 2002) and in 
the current Suzaku XIS observation, which implies a compactness on sizescales of several 
R % . We also know that there is no attenuation of the high energy spectrum below 50keV, 
i.e. £ cut < 511/r 2 < 50keV and hence r < 3. This may also suggest the outflow is more 
compact, i.e. /? w ; ud < 50/?„, than implied for the simple homogeneous case, if the wind is 
launched from e.g. ~ 15i? g . 


Thus we can consider the case where the absorption derives from clumps of matter of 
approximately constant density. In this scenario the maximum distance of the clumps (/? ma x) 
is given by A R/R nmx < 1, where AR is the clump thickness. Now A dI = nAR and n = 
I//£ii^, ax as before; thus we derive the condition that R miiX < L/N h£ and so R tlUiX < 100 7? g 
for A R/R nmx < 1. Thus for R = 20i? gl n ~ 10 9 cm" 3 and A R/R = 0.2. Indeed the fact 
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that both high ionization (log( = 4.9) and lower ionization (partially covering) gas (with 
log £ < 3) is observed in the Suzaku and 2001 XMM-Newton spectra does suggest that the 
absorbing matter is inhomogeneous. The gas may also become less ionized to larger radii, 
especially if the outer layers are shielded from the photoionizing X-ray source. Such matter 
could contribute towards the blue-shifted Lyman a absorption and Civ emission seen in the 
HST/STIS spectrum of PDS 456 (O’Brien et al. 2005), although with columns much smaller 
than those derived here. 


6.3. Outflow Emission and Energetics 

It may be plausible that the outflowing gas produces both the Fe L and K-shell emission. 
The emission could occur via transmission, but given the high column densities involved, 
more likely through reflection off the surface of the wind. To test this, we constructed 
a similar model to the Compton-thick absorber model in Section 5.2, with two absorbing 
layers; one representing a very high ionization outflowing zone and a lower ionization partial 
covering zone (e.g. representing denser clumps within the wind). However in addition to 
the absorption, the model also allows for some reflected emission (via the reflion model) 
from high ionization matter to represent possible scattering off the wind. The illuminating 
continuum is then just a single power-law of T = 2.25 ± 0.03. The reflection component is 
also convolved through a simple Gaussian velocity-broadening function (with A E/E constant 
with energy). This accounts for the line broadening clearly visible in the 7 keV and 1 keV 
emission lines detected by XIS. The model can be expressed as:- 

F{E) = [(CPL x / x PC) + (CPL x (1 - /)) + (REF x GA*^)] x XSTAR x e -^ E > N ».^ (8) 

where CPL is the cut-off powcrlaw with E cxll = 300 keV, PC is the partially covering ab- 
sorber with a covering fraction /, XSTAR is the highly ionized outflowing absorber, REF is 
the reflected emission off partially ionized matter and GA conv corresponds to the Gaussian 
velocity convolution for the reflected emission. 

This model is shown in Figure 9 and the fit parameters are summarised in Table 4 (the 
“outflowing wind” model). A good fit is obtained (* 2 /dof = 217.3/211) and the model well 
reproduces both the absorption in the XIS and HXD spectra as well as the ionized iron L 
and K-shell emission and the general shape of the continuum. The ionization derived for 
the reflector Is high, with !og£ = 3.0 ± 0.3 erg cm s” 1 , which is intermediate between the 
high ionization outflowing absorber and the partial coverer. The ionized reflector represents 
approximately 20% of the broad-band flux of the power-law continuum and hence has R = 
0 . 2 , but contributes less at energies above 10 keV due to the high ionization of this component. 


Thus the scattered component may cover a solid angle corresponding to 6 ~ 0.2 for significant 
emission to be observed in the Suzaku spectrum. We note that a value of b = 0.3 was recently 
derived for the outflow in PG 1211+143, via a P-Cygni like profile to the iron K absorption 
line, suggesting that the kinetic power of the outflow in that quasar is also likely to be similar 
to its total bolometric output (Pounds & Reeves 2009). 

Note that velocity broadening is required for the reflected emission, corresponding to 
35000 km s" 1 , driven by the widths of the emission lines in the spectrum, while some net 
outflow is also preferred (see Table 4). The velocity broadening may arise from the integrated 
emission expected over a wind with appreciable solid angle. A future comparison with 
disk wind models being developed that fully incorporate self-consistent radiative transfer 
calculations for both the emission and the absorption (Sim et al. 2008; Schurch et al. 2009) 
may prove to be interesting. 

As a consistency check we compare the mass outflow rate calculated in Section 6.2 with 
that predicted via transfer of momentum through radiation pressure, assuming that the wind 
is radiatively driven. For such a high ionization outflow, the predominant form of momentum 
transfer to the wind will be through scattering (King &: Pounds 2003), whereby:- 

Mo U ,«W = — (1-0 (9) 

c 

where r ~ 1 as measured in the Suzaku spectrum of PDS456 and Le dd = 2 x 10 47 ergs _1 
(see Section 6.1). Thus from the above M out v out = 4.5 x 10 36 gems" 2 and as v out = 0.26c, 
then the mass outflow rate is M a ut = 6 x 10 26 gs _1 or ~ 10 M© yr~ l . In comparison the mass 
outflow rate calculated in Section 6.2 (equation 5) is A/ 0 ut = 66 x lO^gs" 1 and thus 6 ~ 0.1 
for PDS 456, consistent with the estimate from the above model. 

Nonetheless even if 6 ~ 0.1 then the outflow is energetically significant, with a kinetic 
power of E = 2 x 10 46 ergs -1 . Integrated over a quasar lifetime of 10 8 years and if we assume 
a conservative duty cycle of 10% for the outflow to be active, then the total amount of 
energy released by the wind will be > 5 x 10 6O erg, plausibly exceeding the binding energy of 
a galaxy bulge of mass 10 11 M© and velocity dispersion a — 300 km s" 1 of 5 x 10 ;>a crg. Thus 
such powerful winds could potentially cause significant feedback between the galactic bulge 
and black hole during the quasar phase. 

The duty cycle for the outflow is largely unknown, however we know that the iron 
K absorption was detected in most of the observations of PDS 456 to date (see Figure 7). 
There are also several reported cases of high velocity outflows that are emerging in the 
literature, e.g. Cappi (2006); Reeves et al. (2008) and references therein. The frequency 
of these winds and how the outflow depends on critical parameters (such as accretion rate) 
awaits a more systematic and statistically rigorous survey of such systems in the Suzaku 
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and XMM-Newton archives. Future calorimeter based X-ray spectra with Astro-H and IXO, 
with resolution E/AE > 1000 at 6keV, will hopefully reveal a wealth of information on 
these outflows in the iron K band. 


7. Conclusions 

We have detected the signature of a possible Compton thick wind from PDS 456, a 
high luminosity nearby quasar (z = 0.184) which accretes near to the Eddington limit. The 
Suzaku XIS spectrum shows evidence for statistically significant (at > 99.99% confidence) 
absorption in the Fe K band, which can be modeled by a pair of absorption lines (or one broad 
line) near 9 keV in the quasar rest frame. The large velocity shift, compared to the most 
likely identification of the absorption with FeXXVi (Is - 2 p) at 6.97 keV, implies an outflow 
velocity in the range from 0.26 -0.31c. As has been discussed in detail, an identification with 
other atomic tiansitions is less likely and the large outflow velocity appears to be required. 
The large velocity implies that the mass outflow rate is high, of the order of several tens of 
Solar masses per year, with a corresponding kinetic power of up to ~ 10 47 ergs -1 , close to 
the Eddington-limited luminosity of PDS 456. Such winds could be an important source of 
feedback regulating black hole and bulge growth (King 2003; Di Matteo et al. 2005) during 
the quasar growth phase in the early Universe. 

A tentative hard X-ray detection of PDS 456 above 15keV has also been made (at the 
- 3 <t level) in the HXD/PIN, at a low flux level of 7 x 10" 12 erg cm -2 s" 1 . The hard X-ray 
eriiission requires either a Compton- thick (TVh > 10 24 cm -2 ) partial covering absorber or 
strong reflection, or a combination of both, to account for the apparent hard X-ray excess. 
However at the current low significance level of the detection, a further deep observation of 
PDS 456 with Suzaku is needed to confirm the nature of the hard X-ray excess. 
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Appendix A - Possibility of HXD Contamination 

It is important to discuss the possibility of contaminating X-ray emission from other 
sources within the HXD field of view. Nonetheless there are no known bright (F 2 _ l0 > 
10 erg cm 2 s ! ) contaminating point sources from previous imaging observations of PDS 
456 below 10 keV, e.g. with ASCA GIS, SAX/MECS, XMM-Newton EPIC which have 
similar fields of view to Suzaku/HXD. Above lOkeV, no known hard X-ray sources (other 
than the expected detection of PDS 456 itself) are present within 2° of PDS 456 in the RXTE 
slew survey (Revnivtsev et al. 2004) nor from the Integral/IBIS all sky survey (Krivonos et 
al. 2007a). The nearest known source detected in the Integral all sky survey is the LMXB 
GX9+9, approximately 3° from PDS 456. There is a possibility of some weak Galactic 
diffuse contamination due to the position of PDS 456, at the very edge of the Galactic 
bulge, l = 10.4°, b = +11.2°), due to a population of accreting magnetic CVs. Results from 
the Integral (Krivonos et al. 2007b) or RXTE/PCA (Revnivtsev et al. 2006) Galactic ridge 
surveys show that this contributes < 1 x 10“ 12 erg cm" 2 s" 1 from 15-50 keV, integrated over 
the much smaller HXD/PIN field of view and thus is negligible. 

There is also the contribution of Extragalactic diffuse (AGN) emission from the Cosmic 
hard X-ray Background (CXB), an estimate of which has already been included in the 
background spectrum for the HXD/PIN. The HXD detection of PDS 456 lie's above the level 
expected from the CXB. The flux measured by the HXD, minus the non X-ray background 
(NXB) component (but not the CXB) is 1.52±0,20x 10 11 erg cm" 2 s' 1 (15-50 keV band). In 
this paper we have adopted the expected CXB flux from the HEAO-1 measurement of Gruber 
et al. (1999), which is 7.9 x 10“ 12 ergcnr 2 * * s" 1 from 15-50 keV, thus the net intrinsic flux from 
PDS 456 is 7.3 ± 2.4 x 10" 11 erg cm" 2 s l . The error includes the 1.3% absolute uncertainty in 
the HXD/PIN background, and allows for 10%, uncertainty in the CXB determination. Note 
that the CXB flux measured by Gruber et al. (1999) is also consistent with that obtained by 
Beppo-SAX PDS by Frontera et al. (2007), who obtain a 90% upper bound on the CXB level 
of 6.5 ± 0.2 x 10“ 8 erg cm -2 s~ l sr" 1 (20-50 keV) equivalent to 8.3 ±0.2 x 10" l2 ergcm" 2 s" 1 , 
when normalized to the Suzaku HXD field of view and converted to the 15-50 keV band. We 
note that the Integral IBIS measurement of the CXB, by Churazov et. al. (2007) via Earth 
occultation, is also within 10% of the HEAO-1 measurement. The effect of any anisotropy 
of the CXB should be negligible, fluctuations are thought to be at the level of 5 10% on 
scales of a square degree. As a final check of the CXB level, we also computed the HXD 
flux from the Suzaku observation of the Lockman Hole (LH), performed on 3 6 May 2007. 
The data were reduced as described in Section 2, resulting in 84.3 ks of net HXD exposure. 
The HXD/PIN spectrum of PDS 456 lies significantly above that of the LH emission, while 
the LH flux is consistent with the expected CXB level (0.9 i 0.2 x 10" 11 erg cm" 2 s" 1 ). Thus 
uncertainty in the CXB subtraction does not significantly effect the detection of PDS 456. 
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One final possibility is that a bright, serendipitous AGN is located within the HXD 
field of view. From the the Swift BAT hard X-ray selected AGN logN - logS (Tueller et 
ah 2008), we would randomly expect < 5 x 10~ 3 AGN of flux ~ 1 x 10~ 11 erg cm” 2 s -1 
per HXD field. The actual probability Is likely much lower, as there are no other bright 
AGN candidates in the soft X-ray (< lOkeV) or optical fields of view. Thus in order to 
contaminate the HXD/ PIN field of view, a background AGN would itself have to be heavily 
absorbed (Nu > 10 21 cm" 2 ) below lOkeV (so that F 2 ~ io <; 10“ 13 ergcm _2 s _1 , as there are 
no known bright source in the PDS 456 field below 10 keV) and yet its spectrum must 
rise by a factor of xlOO above 10 keV. Statistically speaking, < 5% of the BAT selected 
AGN have Nu > 10 21 cm 2 , meaning that the coincidence probability is ~ 10~ 4 . The only 
other possibility is of a flaring object, such as a blazar, brightening during the observation, 
although the HXD lightcurve is flat (Figure 1), which appears to exclude this and also the 
object would have to be located outside the imaging field of view of the XIS. 

Furthermore, PDS 456 was also detected previously in both 1998 (Reeves et al. 2000) and 
2001 with RXTE, with a hard X-ray flux of F 3 _ 15 = 0.7 — 0.9 x 10~ u ergcm~ 2 s _1 . PDS456 
was also more recently detected in a 35 ks worth of RXTE pointings in 2008 (see Table 1 for 
a list of observations) . PDS 456 was not detected in a previous observation with Bcppo-SAX 
PDS (Feb 2001, 71 ks PDS exposure), however the upper limit obtained is consistent with 
the above HXD analysis (F 2 o- u>o < 1*5 x 10~ u erg cm" 2 s -1 ). Indeed although the RXTE 
observations and the Suzakju HXD observation are not simultaneous (see Table 1), the time- 
averaged RXTE/PCA spectrum extrapolates very well to the HXD/PIN measurement (see 
Figure 10). Therefore it appears likely that the detected hard X-ray emission is due to the 
quasar and not; due to any; contaminating sources in the field of view. 
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Fig. 1. — (a) Comparison between the background A model and the Earth occult at ion data 
for HXD/PIN (from 15— 50keV, where the light curve represents the ratio between the (Earth 
- Bgd A) / Bgd A rate. The background A model underpredicts the Earth background by 
8.8 ± 1.0%. (b) Comparison between the background D model and the Earth occupation 
data, where the light curve represents the ratio between the (Earth — Bgd D) / Bgd D rate. 
The background D model slightly underpredicts the Earth background by 2.4±1.0% (c) The 
net 15-50 keV HXD/PIN lightcurve using the background D model, the CXB background 
level (0.17 ctss -1 ) is marked by a horizontal line. 
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Fig. 2. Upper panel - Suzaku XIS spectra of PDS 456, plotted as a ratio to a simple 
absorbed power-law (r = 2.2). Note XISl is shown in red, XIS FI in black. Lower panel 
- rest frame iron K-shell absorption in PDS 456 (against the best fit, broken power-law 
continuum). Suzaku XIS FI is shown in black, XMM-Newton EPIC-pn (2001) in red. The 
green line represents the expected position of the Fexxvi Is - 2 p line (6.97 keV), while the 
blue lines represent the centroids of the defected absorption lines in XIS. 
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Observed Energy (keV) 

Fig. 3. Suzaku XIS FI PDS 456 spectrum versus the background spectrum in the iron 
K band. The source spectrum prior to background subtraction in shown in red, the net 
source spectrum after background subtraction is in black and the background spectrum is 
in green. The black line shows the broken powerlaw continuum model, folded through the 
XIS response. It can be seen that the background subtraction has little effect in the iron K 
band and on the absorption features at 7.7 keV and 8.2 keV (observed frame). Overall the 
level of the XIS background is low and the only strong background line arises from Ni Ka 
at 7.45 keV, which is not coincident with the centroids of the absorption lines from PDS 456. 
Note the normalization of the background has been rescaled to that of the XIS source region 


size. 
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Fig. 4. Comparison between outflowing Xstar models fitted to the Suzaku XIS-FI, showing 
the iron K band. Panel (a) shows a fit with one Xstar- zone, outflowing with a single velocity 
of v ou t = —0.30c and assuming a turbulence velocity of t; turb = 10000 km s -1 (model A in 
Section 4.2). The absorption model is shown by the red line, the 8keV absorption (observed 
frame) is modeled by a blend of Fexxv and Fexxv Is — 2 p resonance lines, while iron 
K emission is shown in green. Panel (b) shows a similar model as per (a), but with two 
outflowing zones of v mi = f- 0.26c and —0.31c (model B) Panel (c) shows a similar model as 
per panel (a), but with Gaussian velocity broadening included self-consistently for the iron 
I< and L- shell emission and absorption (model C). The blue line shows the model without 
broadening. Further details are given in Section 4.2 and parameters for models A, B and C 
are listed in Table 3. 



Fig. 5. Xstar fit to the Suzaku XIS FI data as per Figure 4, but with no net outflow 
velocity and assuming a very low turbulence velocity (100 kms" 1 ); e.g. model D in Section 
4.2 and Table 3. Here the absorption at 8 keV in the observed frame is modeled by bound- 
free edges from Fexxv and Fexxvi, from a high column density absorber (i.e. r = 1), while 
narrow Is — 2 p resonance lines from Fexxv and Fexxvi appear near 6keV in the model 
(but not in the data). The fit is statistically worse than for the outflowing models shown in 
Figure 4 and such a small turbulence may be physically implausible. 
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Fig. 6. The Suzaku XIS and HXD spectra of PDS 456, plotted in i ;F V flux units. A 
factor of x3 excess is seen in the HXD, which may suggest that a Compton thick absorber, 
with Nu > 10 24 cm -2 exists towards this optical type I quasar. The red line shows the 
extrapolation of the best-fit broken-powerlaw continuum extrapolated to higher energies. 
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Fig. 7. Suzaku XIS 4- HXD spectrum of PDS 456, fitted with an ionized disk reflection 
model, as described in Section 5.1. Data are shown in black, the blue line shows the blurred 
reflection component, green shows the intrinsic power-law continuum and red the total emis- 
sion. Note the HXD/PIN spectrum is more finely binned than in Figure 6, at ~ 2a per bin 
above the background level. 
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Fig- 8. The current observations of PDS 456, shown as vF v spectra, where we have plotted 
the data compared to a ratio with an absorbed (N G ai = 2x 10 21 cm -2 ) power-law of T = 2 
at a mean 2 lOkeV flux level of 5 x 10 12 ergcm~ 2 s“ l . Shown are RXTE PCA (black, 
1998), XMM-Ncwton (2001, green), Suzaku XIS (2007, red), Chandra/HETG (2003, blue) 
and ASCA (1998, orange dotted). Strong variability is seen below 10 keV, from RXTE and 
XMM at high fluxes, down to Chandra and ASCA as the lowest /hardest observations. Note 
the 8 keV (observed) Fc K-shell absorption in the XMM 2001, Suzaku-2007 and RXTE 1998 
observations. A summary of all these observations is shown in Table 1. 



Fig. 9.— Model fitted to the Suzaku spectrum, representative of absorption and emission 
from a high column density disk wind, as described in Section 6.3, with model parameters 
listed in Table 4. Datapoints are shown in black, the dotted green line represents the power- 
law (T = 2.2) continuum that is absorbed only by the very high ionization fast outflowing 
zone, the dark blue line is the continuum obscured by the partial covering absorber which 
covers 80% of the X-ray source, the light-blue line is the reflected emission off the Compton- 
thick partially ionized wind and red is the total emission from all components. 
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Fig. 10. A comparison between the time averaged RXTE/PCA spectrum between 3 - 
18keV and the HXD/PIN Suzaku spectrum from 15 50keV. The plot shows the comparison 
of the data to a T = 2.2 continuum with Galactic absorption, with no re-normalization of 
either spectrum, while the lower panel shows the residuals between the data and the model. 
Although the RXTE/PCA data is non simultaneous with the HXD (see Table 1 for a list 
of observations), the PCA spectrum below 18keV extrapolates well to the HXD/PIN data 
above 15 keV. Also notice the deep absorption trough centered at 8 keV in the RXTE/PCA 
spectrum, consistent with the absorption in the 2007 Suzaku XIS and 2001 XMM-Newton 
data. 


Table 1. Summary of PDS 456 Observations 


Telescope 

Mode" 

Start Date/Tirae** 

End Date/Time 6 

Exposure* 

F 2 - w d 

Suzaku XIS 

XIS nom 

2007/02/24 17:58:04 

2007/03/01 00:51:14 

179.6 

3.7 

Suzaku HXD/PIN 




164.8 


XMM-Newton pn 

FW/med 

2001/02/26 10:46:57 

2001/02/26 22:00:41 

36.4 

6.9 

Chandra/HETG 


2003/05/07 03:29:56 

2003/05/08 20:08:33 

142.8 

3.8 

ASCA/SIS 


1998/03/07 09:06:36 

1998/03/08 14:58:42 

42.8 

4.0 

RXTE/PCA 


1998/03/07 17:07:44 

1998/03/10 07:50:56 

92.1 

8.1 



2001/02/23 04:20:32 

2001/03/11 05:59:44 

168.3 

10.4 



2008/02/22 21:00:32 

2008/11/14 08:17:36 

34.7 

69 

Beppo-SAX/PDS 


2001/02/26 09:50:03 

2001/03/03 08:57:30 

71.0 



instrument mode or pointing. XIS nom = XIS nominal pointing; FW = Full Window; mod = 
medium filter. 

b Observation Start/End times tire in UT. 

c Net exposure time, after screening and deadtime correction, in ks. 
d Absorbed flux in the 2-10 keV band, units x 10“ 12 erg cm 2 s' 1 
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Table 2. Spectral Parameters to Suzaku XIS 


Model Component 

Fit Parameter 

Value 

& x 2» 

1. Broken Power-law Continuum 6 

rWt 

2.44 ± 0.04 



Mm*ak 

2.0 ± 0.2 



r hard 

2.24 ± 0.02 



Nr PL 

(2.4 ±0.1) x 10 ~ 3 



Fo.r>- mfccV 

7.2 x 10 “ 12 


2 . Galactic Absorption 6 ' 

Mt 

( 2.2 ± 0 . 2 ) x 10 21 


3. Soft X-ray Emission Lines'* 

Mine 

0.91 ±0.01 

174.4 


Line Flux 

(2.0 ± 0.4) x lO " 4 



EW 

42 ± 10 



a 

0.075 ± 0.015 


Second line 

M ine 

1.16 ± 0.02 

17.7 


Line Flux 

(1.4 ±0.6) x 10 " 5 



EW 

6 ± 3 



a 

< 0.05 


4. Iron K Absorption Lines'* 

Mine 

9.09 ± 0.05 

41.1 


Line Flux 

—3.0±J;f x 10 ~ 6 



EW 

137 ± 59 


Second line 

Mine 

9.64 ± 0.08 

29.0 


Line Flux 

-(2.5 ± 1.0) x 10 “ 6 



EW 

135 ± 55 



a 

0. 105^0^1 ( tie d to fi rst ^ ine ) 


5. Iron K Emission Line'* 

Mine 

7.1 ±0.2 

27.7 


Line Flux 

(5.1 ±2.0) x 10 “ 6 



EW 

132 ± 52 eV 



a 

o-Kijs 


Fit statistics ' 7 

X 2 /dof 

431.5/374 



Null probability 

2.1 x 10" 2 



Table 3. Xstar Photoionization Model Fits to Suzaku XIS 


Model® 

jV 

log£ c 

aturb^ 

Vout 6 

Pbroad^ 

X 2 /dof 3 

Null Prob 5 

A 

0 99 +°- 11 

4.1 ±0.2 

10000 

-0.30 ± 0.02c 


224.9/204 

0.151 

B 

1.3 ± 0.6 

4 q+0-9 

4y -0.5 

10000 

-0.26 ± 0 . 02 c 


207.3/203 

0.394 


tied 

tied 

tied 

-0.31 ± 0.02c 




C 

1.5 ± 0-6 

4.4 ± 0.4 

10000 

-0.29 ± 0.02c 

26500 ± 7000 

227.9/205 

0.131 

D 

9 0 + 1.2 
Z - 6 -Q.7 


100 

0 


254.0/208 

1.6 x 10“ 2 


a Models A, B, C and D as defined in the text (section 4.2). Model A contains a single pho- 
toionized absorber with one outflow velocity; model B contains two absorbers with two seperate 
outflow velocities; model C contains one outflow velocity with equal kinematical broadening in emis- 
sion/absorption; model D has zero outflow velocity and low (lOOkms -1 ) turbulence. 

b Column density, in units xlO 24 cm -2 . 
c Ionization parameter, units ergs cm s” 1 . 

^Turbulence velocity width (la), in units kms -1 . 
e Outflow velocity, in units of c. 

* Gaussian velocity broadening (FWHM), in units kms” 1 . 

g Reduced chi-squared (y 2 /dof) and null hypothesis probability for spectral fit. 


‘‘Improvement in \ 2 to fit after adding model component. 

b Broken power-law continuum parameters. Thard and r w) f t , soft and hard photon indi- 
cies; inbreak i break energy in units keV; Mjpl> broken powerlaw normalization in units pho- 
tons cm 2 s 1 at IkeV; Fq lofccvs absorbed continuum flux in units ergs cm -2 s -1 . 

c Local Galactic absorption (at z = 0 ), units cm -2 . 

d Emission/ Absorption line parameters. Minet line energy in units keV; line flux in units pho- 
tons cm 2 s equivalent width (EW) in units eV; a is la line width in units keV. Parameters 

given in quasar rest- frame. 

‘ Reduced chi-squared ( \' 2 / dof) and null hypothesis probability for spectral fit. 
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Table 4. Outflowing Wind Model Parameters 


Model Component 

Fit Parameter 

Value 

1. Power- law Continuum 

r 

2.25 ±0.05 


•^abs' 1 

(8.1 ±2.3) x 10~ 3 


JV„„ abs a 

(2.0 ±0.1) x 10~ 3 


n /; 

r 0.5- 10 keV 

7.2 x 10" 12 


Fio-SOfceV^ 

7.8 x 10- 12 

2. Outflowing Fe K-shell Absorber 

N U C 

(1.3 ±0.5) x 10 21 


log^ 

4 q+o 9 
4 - y -o.5 


iwC 

-0.26 ± 0.02c 


t,' ou t2 e 

-0.31 ± 0.02c 

3. Compton-thick Hard X-ray Absorber 

N h c 

> 2.5 x 10 24 


log^ 

2.5 ±0.4 


f CUV f 

0.80^ 7 

4. Ionized Reflection 

\ogi d 

3.0 ± 0.3 


Vbroad^ 

35000 ± 17000 


’’refi' 1 

-42600 ± 5000 


R 

0.2 

Fit statistics 

X 2 /doP 

217.3/211 


Null probability 5 

0.367 


a Powerlaw normalization in units photons on -2 s -1 at IkeV; Fo.S-iot-cV'- N^ s 
refers to the power-law component that is absorbed by the hard X-ray absorber, 
N unabs is the component not covered by the hard X-ray absorber. 

b Measured flux in units ergs cm” 2 s -1 . 
c Column density, units cm" 2 . 
d Ionization parameter, units ergs cm s -1 . 

'“Iron K absorber outflow velocity, in units of c. 
f Absorber covering fraction 

- FWHM Gaussian velocity broadening in units kins -1 . 
h Net outflow velocity of reflector in units kms -1 . 

'Ratio of reflected emission to power- law continuum. 

■■Reduced chi-squared (y 2 /dof) and null hypothesis probability for spectral fit. 


